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1 .  SUMMARY 

A  time-dependent  Monte  Carlo  computer  code  which  simulates  the  far-field 
interaction  between  a  rocket  exhaust  plume  and  an  ambient  flowing  gas  was 
developed.  Computer  calculations  were  carried  out  for  a  rocket  motor  firing 
in  both  continuous  and  pulsed  modes.  The  pulsed  mode  calculations  employed 
the  Baum,  Kolb,  Rieger,  and  Sutton  Radiation  Model ^  for  the  atomic  physics 
of  the  H^O  molecule,  which  provides  estimates  of  radiant  flux  from  the  inter¬ 
action  region.  Fran  the  pulsed  mode  calculations,  it  was  found  that,  for  the 
assumed  rocket  motor  properties  and  operating  characteristics  of  600  pounds 
thrust  and  .007  seconds  bum  time: 

1.  At  altitudes  above  200  km,  for  times  less  than  3  exhaust  molecule 
mean  free  times,  and  for  angles  between  relative  wind  and  sensor 
look  axis  of  90°  or  less: 

a)  the  flow  is  in  the  "free  streaming  shell  regime" 

b)  the  order  of  magnitude  of  the  power  on  a  detector  can  be 
estimated  by  a  simple  calculation  --  the  dominant  contri¬ 
bution  comes  directly  from  the  free  molecular  interaction 
between  the  shell  of  exhaust  gas  and  the  ambient  gas 

c)  the  power  on  a  detector  can  be  scaled  simply  with  magni¬ 
tude  and  orientation  of  the  relative  wind  by  taking  the 
vector  sum  of  U.  and  Uexha||st. 

2.  The  time  for  the  disturbance  to  be  blown  away  is  sensitive  to  the 
orientation  of  the  relative  wind. 


H.  Baum,  C.  Kolb,  T.  Rieger,  E.  Sutton,  A.  Germeles;  "(S)  Mater  Vapor 
Rotational  Band  Radiation  From  High  Altitude  Rocket  Plumes  (U) "  Aerodyne 
Research,  Inc.  Report  No.  RN-6,  March  1972. 
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2.  INTRODUCTION  AND  BACKGROUND 

At  high  altitudes  in  the  mixing  zone  between  the  highly  expanded  portion 
of  a  rocket  exhaust  plume  (the  plume  "far  field")  and  an  ambient  supersonic 
stream,  the  flow  is  strongly  influenced  by  non-equilibrium  molecular  transport 
processes.  Chemical  reactions  occurring  in  the  flow  may  not  be  describable 
by  thermal  rate  coefficients  and  proper  description  of  radiation  from  the 
mixing  zone  may  require  an  accurate  determination  of  the  molecular  velocity 
distribution. 

A  method  which  has  been  demonstrated  to  be  very  successful  in  describing 
such  flows  is  the  Monte  Carlo  Simulation  Technique,  in  which  the  molecular 
motion  is  simulated  on  a  digital  computer  using  the  Monte  Carlo  random  sampling 
method.  The  technique  is  well  suited  to  treat  geometrically  complex,  time 
dependent,  three-dimensional  flows.  Recent  extensions  to  the  technique  enable 
chemical  reactions  and  radiation  to  be  described  in  the  simulation. 

A  Monte  Carlo  code  has  been  developed  to  describe  the  radiation  due  to 
the  time-dependent,  three-dimensional  interaction  in  the  far  field  between 
an  ambient  flowing  gas  and  a  cloud  of  exhaust  gas  molecules  emitted  by  a 
rocket  firing  pulse  of  specified  duration.  The  computer  code  capabilities 
include  the  following: 

a)  The  angle  between  the  rocket  exhaust  axis  and  the  relative  wind  is 
an  input  variable,  to  enable  the  effect  of  orientation  to  be 
identified. 

b)  As  initial  condition  for  the  interaction  calculation,  the  cloud  of 
exhaust  gas  molecules  is  assumed  to  be  "inertia  dominated"  and 
hence  is  prescribed  to  form  a  source-like  collisionless  flow  which 
can  be  well  approximated  using  inviscid  method-of-characteri sties 
theory. 

c)  A  "general  radiation  model"  is  incorporated  In  the  code.  The  reac¬ 
tive  collision  cross  sections  to  create  and  destroy  radiatively- 
active  species  are  assumed  as  input  data, as  a  function  of  impact 
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parameter  and  relative  kinetic  energy.  The  mean  radiative  life¬ 
time  is  assumed.  Using  this  information,  during  the  calculation 
of  the  flov/field,  collisional  excitation  of  molecules  is  allowed 
to  occur  with  the  appropriate  probability.  For  each  newly  excited 
molecule,  a  radiative  lifetime  is  sampled  from  the  appropriate 
distribution.  Each  excited  molecule  is  tracked  in  the  flow, 
checking  each  subsequent  collision  for  the  possibility  of  deactiv¬ 
ation,  until  the  time  to  radiate  is  reached.  The  resulting  photon 
flux,  as  a  function  of  location  in  space,  is  recorded  during  the 
calculation. 

d)  As  an  alternative  to  the  radiation  model  described  in  c)  above, 

the  "effective  radiation  rate  constant"  from  the  Baum,  Kolb,  Reiger, 
and  Sutton  Radiation  Model ^  for  ^0  was  also  incorporated  into  the 
computer  program.  This  "rate  constant"  is  essentially  the  average 
energy  released  per  collision  and  is  given  as  *  function  of  relative 
velocity.  Counters  for  each  small  region  of  pnysical  space  were  set 
up  in  the  code  and  were  incremented  at  each  collision  by  the  appro¬ 
priate  amount.  The  emitted  power  from  each  region  is  then  found  by 
differentiating  the  cumulative  energy  curve. 


21 222-6001 -R0-00 
Page  4 


3.  DISCUSSION  OF  RESULTS 


3.1  General 

The  time-dependent  Monte  Carlo  computer  code  to  simulate  far-field 
plume  interaction  was  completed  at  TRW  in  October  1971.  Preliminary  cal¬ 
culations  to  verify  the  code  were  carried  out  in  October  at  the  Safeguard 
computing  facility  in  Huntsville,  Alabama.  Subsequently,  a  set  of  produc¬ 
tion  calculations  was  carried  out  at  Huntsville.  The  conditions  for  the 
calculations  are  shown  in  Figure  1. 

Typical  rocket  motor  characteristics  were  taken  to  be  600  pounds  thrust 
and  7  mill i-seconds  burn  time.  The  vacuum  expansion  flowfield  used  for 
initial  conditions  at  the  inner  boundary  of  the  field  of  view  was  generated 
by  AVCO  Everett  Research  Laboratory  using  the  method  of  F.  P.  Boynton. 

3.2  Preliminary  Steady  Flow  Calculations 

Two  steady  flow  calculations,  at  simulated  altitudes  of  140  km  and 
200  km,  were  carried  out.  These  were  intended  to  investigate  possible 
differences  between  the  far-field  interaction  due  to  a  rocket  motor  firing 
in  a  continuous  mode  and  one  firing  in  a  pulsed  mode. 

Stagnation  line  profiles  of  the  relative  fraction  of  exhaust  molecules 
and  ambient  molecules  at  200  km  are  shown  in  Figure  2.  The  parameters  char¬ 
acteristic  of  our  assumed  maneuvering  rocket  motor  were  used  for  this  cal¬ 
culation,  except  that  the  bum  time  was  made  infinite,  and  the  relative  wind 
was  aligned  anti-parallel  to  the  thrust  axis;  that  is,  the  rocket  motor 
is  in  retro-fire.  The  fractions  of  exhaust  molecules  and  ambient  molecules 
which  have  been  scattered  from  the  undisturbed  exhaust  or  ambient  state 
either  by  collisions  with  other  scattered  molecules  or  exhaust-ambient 
collisions  are  also  shown  in  Figure  2.  It  can  be  seen  that  the  shock  layer 
is  well  "merged";  that  is,  external  and  internal  shock  waves  are  not  apparent 
as  distinct  regions. 


MONTE  CARLO  CALCULATIONS  CARRIED  OUT 
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FIGURE  1 


Stagnation  Line 
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Comparison  of  these  steady  flow  results  and  the  unsteady  flow  results 
to  be  discussed  below  revealed  substantial  quantitative  differences.  Steady 
flow  results  thus  would  not  be  useful  to  make  predictions  for  a  motor  firing 
in  the  pulsed  mode. 

3.3  Unsteady  Flow  Calculations 

A  tabulation  of  some  characteristic  flow  interaction  parameters  for  the 
assumed  rocket  motor  firing  intermittently  is  shown  in  Figure  3.  Values  of 
the  ratio  of  shell  width  to  exhaust  molecule  mean  free  path,  &£/Xef,  are  shown 
there.  It  seems  likely  that  for  altitudes  of  200  km  and  above,  where  Ai/x^ 
is  «  1,  the  shell  of  gas  emitted  by  the  rocket  is  too  thin  for  there  to  be 
any  significant  "cascading"  of  collisions  within  the  shell  of  exhaust  gas. 

One  would  then  expect  the  gas  density  in  the  shell  to  simply  attenuate  with 
distance  from  the  nozzle  due  to  geometry  and  due  to  collisions  with  ambient 
molecules.  For  altitudes  between  140  km  and  200  km  a  conclusion  based  upon 
such  a  simple  analysis  must  be  more  tentative;  however,  it  seems  likely  that 
somewhere  between  140  km  and  200  km  the  effect  of  "cascading"  will  become 
appreciable.  It  is  not  possible  to  make  a  similar  a  priori  estimate  for  the 
interaction  with  the  atmosphere  of  the  molecules  scattered  from  the  shell. 

3.3.1  Unsteady  Calculation  at  90°  Angle  of  Attack 

An  unsteady  calculation  at  200  km  altitude  with  the  relative  wind 
aligned  anti-parallel  to  the  thrust  axis  (90°  angle  of  attack)  was  carried 
out  using  the  "general  radiation  model"  described  in  2c).  In  Figure  4,  the 
relative  composition  of  the  gas  along  the  thrust  axis  is  shown  at  60  milli¬ 
seconds,  or  one  mean  free  time,  after  the  start  of  the  calculation.  Five 
species  are  identified,  including  exhaust  molecules  which  have  been  excited 
to  an  IR  active  state.  The  gas  behind  the  shell  is  seen  to  be  composed  of 
approximately  70  percent  undisturbed  free  stream  molecules  and  25  percent 
scattered  free  stream  and  exhaust  molecules.  Excited  molecules  form  about 
5  percent  of  the  total.  It  should  be  pointed  out  that  the  "general  radiation 
model"  used  in  the  code  is  an  extremely  simple  one,  utilizing  a  single  step 
excitation.  A  more  realistic  radiation  model  was  supplied  by  ARI^  and  used 
in  a  subsequent  calculation,  to  be  discussed  in  3.3.2. 


CHARACTERISTIC  PARAMETERS 
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FIGURE  3 
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Fic.re  5  shows  the  cumulative  radiative  flux  in  two  locations  in  the 
fie./  fi  I d,  a 1  ong  the  thrust  axis  and  look  axis,  at  60  ms  and  110  ms  after 
the  st  (rt  o'  the  calculation.  The  ordinate  on  the  figure  is  the  number  of 
phot  ns  produced  per  local  free  stream  molec  <e.  Due  to  the  simplicity  of 
thr  radiation  model  and  uncertainty  in  the  '  gritude  of  the  excitation 
cu  ss  section,  it  is  not  known  how  much  qua  litative  significance  can  be 
attached  to  the  radiation  information  presented  in  Figure  5.  Probably  the 
most  significant  feature  is  that  the  ordinate  value  never  exceeds  one. 

Since  the  excitation  cross  section  used  in  the  calculation  was  equal  to 
gas  kinetic  (reduced  by  a  geometric  factor),  this  means  that  each  free 
stream  molecule  excites  at  most  one  exhaust  molecule.  Hence,  it  is  clear 
that  there  is  no  "cascading"  occurring  in  the  shell  of  gas;  molecules  are 
simply  knocked  out  of  the  s,  >11  to  form  a  trail  of  scattered  species  behind 
the  shell . 

3.3.2  Unsteady  Calculation  at  0°  Angle  of  Attack 

An  unsteady  calculation  at  200  km  altitude  with  the  relative  wind 
aligned  anti-parallel  to  the  sensor  look  axis  (0°  angle  of  attack)  was  also 
carried  out.  The  flow  field  for  this  case  was  found  to  be  similar  to  that 
found  at  90°  angle  of  attack,  that  is,  the  shell  of  emitted  exhaust  gas  is 
simply  attenuated  due  to  geometry  and  due  to  a  free  molecular  interaction 
with  the  atmosphere,  leaving  a  trail  of  scattered  molecules  behind  it. 

This  calculation  was  terminated  at  3  exhaust  molecule  mean  free  times 
(«180  ms)  after  the  leading  edge  of  the  shell  crossed  the  inner  boundary  of 
the  physical  field  of  view.  (The  exhaust  molecule  mean  free  time  = 

+  vex)).  This  crossing  would  occur  approximately  30  ms  after 
ignition.  The  calculation  at  90°  angle  of  attack  described  In  3.3.1  was 
carried  out  to  5  mean  free  times.  In  both  cases,  scattered  exhaust  molecules 
were  found  to  exist  In  significant  quantity  along  the  sensor  look  axis  at  the 
and  of  the  calculation,  to  a  greater  extent  at  90°  angle  of  attack  than  at  0° 
Thus  the  time  for  the  entire  disturbance  to  be  completely  blown  away  by  the 
ambient  wind  Is  certainly  greater  than  3  mean  free  times  for  all  relative 
wind  orientations  and  will  vary  with  the  orientation  angle.  One  can 


90°  Angle  of  Attack 


meters  from  nozzle 
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conservatively  say  that  for  times  less  than  3  mean  free  times,  the  flow  is 
in  the  "free  streaming  shell"  regime,  as  shown  in  Figure  6. 

To  obtain  a  realistic  estimate  of  the  radiance  from  the  plume  and  the 
power  on  the  detector  the  calculation  described  above  also  included  the 
radiation  model  for  HgO  from  Reference  1.  A  curve  of  the  variation  with  energy 
of  encounter  of  the  average  energy  emitted  per  collision  for  this  model  is 
shown  in  Figure  7.  This  information  was  supplied  to  TRW  by  ARI.  An  analytic 
fit  to  the  ARI  curve  was  incorporated  in  the  Monte  Carlo  computer  program, 
and  at  each  collision  a  counter  was  incremented  by  the  energy  released  appro¬ 
priate  to  the  relative  velocity  for  the  collision,  as  determined  from  the 
ARI  curve.  The  physical  space  surrounding  the  rocket  was  divided  into  small 
incremental  volumes  (called  cells)  each  with  its  own  emitted -energy  counter. 
Typical  time  histories  of  the  cumulative  emitted  energy  are  shown  in  Figure  8 
for  two  incremental  volumes  located  at  r  =  35  meters  and  r  =  45  meters, 
where  r  is  the  distance  from  the  rocket  motor  along  the  sensor  look  axis. 

To  obtain  the  emitted  power  from  each  volume  the  derivative  is  taken.  The 
rapid  initial  rise  for  the  two  histories  in  Figure  8  is  due  to  the  passage 
of  the  shell  of  unattenuated  exhaust  gas.  The  contributions  to  the  cumulative 
emitted  energy  are  extremely  small  thereafter.  Thus,  the  dominant  contri¬ 
bution  to  the  radiant  power  emitted  per  unit  volume  is  clearly  due  to  the 
free  streaming  shell  interaction.  The  total  power  on  the  sensor  is  obtained 
by  the  linear  integration  with  distance  of  the  power  per  unit  volume  along 
the  sensor  look  axis.  For  time  less  than  3  mean  times,  the  contribution 
of  the  "hot"  gas  behind  the  shell,  which  is  comprised  of  the  molecules 
scattered  from  the  shell,  was  found  to  be  negligible.  The  Monte  Carlo  cal¬ 
culation  of  power  on  the  sensor  is  shown  in  Figure  9  and  is  seen  to  be  in 
close  agreement  with  a  free  molecular  hand  calculation,  using  the  Baum,  etal 
model.  The  exhaust  gas  density  for  this  hand  calculation  was  taken  to  be 
the  vacuum  expansion  density  reduced  exponentially  with  distance  to  account 
for  collisions. 

Because  of  the  free  molecular  nature  of  the  interaction,  it  is  possible 
to  scale  the  results  in  Figure  9  to  other  altitudes  and  flight  velocities. 
Normalizing  power  by  the  cube  of  the  ambient  density,  normalizing  time  by 
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FIGURE  7 
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the  exhaust  molecule  mean  free  time,  and  adjusting  the  power  for  other 
velocities  by  an  integration  over  the  free  molecule  relative  velocity  dis¬ 
tribution  using  the  analytic  fit  to  the  Baum,  etal^  curve  in  Figure  7,  a 
generalized  curve  of  power  on  the  sensor  which  is  valid  within  the  "thin 
shell  streaming  regime"  is  shown  in  Figure  10. 
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FIGURE  10 
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4.  CONCLUSIONS 

For  rocket  motor  properties  and  operating  characteristics  of  600 
pounds  thrust  and  7  milli-seconds  burn  time,  it  was  found  that: 

1.  At  200  km  altitude,  for  time  -  3  mean  free  times  (-  .2  secs), 
and  for  angles  between  relative  wind  and  sensor  look  axis  of 
90°  or  less: 

a)  the  flow  is  in  the  "free  streaming  shell  regime" 

b)  the  order  of  magnitude  of  the  power  on  the  detector  can 
be  estimated  by  a  simple  calculation  —  the  dominant  con¬ 
tribution  comes  directly  from  a  free  molecular  inter¬ 
action  between  the  shell  of  exhaust  gas  and  the  ambient 
gas 

c)  the  power  on  a  detector  can  be  scaled  simply  with  magni¬ 
tude  and  orientation  of  relative  wind  by  taking  the  vector 
sum  of  and  Uex^au$t  and  then  employing  the  effective 
radiation  rate  constant  from  Reference  1. 

2.  Conclusions  in  1)  above  apply  to  higher  altitudes  with  simple 
scaling. 

3.  The  time  for  the  disturbance  to  be  blown  away  is  sensitive  to  the 
orientation  of  the  relative  wind,  but  is  greater  than  3  mean  free 
times  for  all  orientations. 


